Quantum Phase Transition in Quantum Dot Trimers 
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We investigate a system of three tunnel-coupled semiconductor quantum dots in a triangular ge- 
ometry, one of which is connected to a metallic lead, in the regime where each dot is essentially 
singly occupied. Both ferro- and antiferromagnetic spin-^ Kondo regimes, separated by a quantum 
phase transition, are shown to arise on tuning the interdot tunnel couplings and should be acces- 
sible experimentally. Even in the ferromagnetically-coupled local moment phase, the Kondo effect 
emerges in the vicinity of the transition at finite temperatures. Physical arguments and numerical 
renormalization group techniques are used to obtain a detailed understanding of the problem. 

PACS numbers: 71.27.+a, 72.15.Qm, 73.63.Kv 



I. INTRODUCTION 

Nanofabrication techniques developed over the last 
decade 1-8 , together with atomic-scale manipulation us- 
ing scanning tunneling microscopy 9,10 , have sparked in- 
tense interest in novel mesoscopic devices where strong 
electron correlation and many-body effects play a central 
role 11 . The classic spin- ^ Kondo effect 12 - in which a sin- 
gle spin is screened by antiferromagnetic coupling to con- 
duction electrons in an attached metallic lead - has been 
observed in odd-electron semiconductor quantum dots, 
single molecule dots, and adatoms on metallic surfaces 11 ; 
although its ferromagnetic counterpart, where the spin 
remains asymptotically free, has yet to be reported exper- 
imentally. In coupled dot devices - the 'molecular' ana- 
logue of single quantum dots viewed as artificial atoms 2 
- exquisite experimental control is now available over ge- 
ometry, capacitance, and tunnel- couplings of the dots 5,6 . 
Both spin and internal, orbital degrees of freedom - and 
hence the interplay between the two - are important in 
such systems. This leads to greater diversity in poten- 
tially observable correlated electron behaviour on cou- 
pling to metallic leads, as evident from a wide range of 
theoretical studies of double (e.g. Refs. 13-21) and triple 
(e.g. Refs. 22-29) quantum dot systems. 

Motivated in part by recent experiments involving 
triple dot devices 4-7 , we consider here a system of three, 
mutually tunnel-coupled single- level quantum dots, one 
of which is connected to a metallic lead: a triple quantum 
dot (TQD) ring structure, the simplest to exhibit frus- 
tration. We focus on the TQD in the 3-electron Coulomb 
blockade valley, and study its evolution as a function of 
the interdot tunnel couplings, using both pcrturbative ar- 
guments and the full density matrix 30,31 formulation of 
Wilson's numerical renormalization group (NRG) tech- 
nique 32,33 (for a recent review, see Ref. 34). A rich range 
of behaviour is found to occur. Both antiferromagnetic 
and ferromagnetic Kondo physics are shown to arise in 
the system - with the two distinct ground states sep- 
arated by a quantum phase transition - and should be 
experimentally accessible via side-gate control of the tun- 
nel couplings. The zero-bias differential conductance (G) 
through the dots is shown to drop discontinuously across 



the transition at zero temperature, from the unitarity 
limit of G/Gq = 1 in the strong coupling antiferromag- 
netic phase (with Go = 2e 2 /h the conductance quan- 
tum) to G/Gq ~ in the weak coupling, local moment 
phase. However in a certain temperature window in the 
vicinity of the transition, the conductance is found to 
be controlled by the transition fixed point separating the 
two ground state phases, comprising both a Kondo sin- 
glet state and a residual local moment; in particular such 
that increasing temperature in the local moment phase 
actually revives the antiferromagnetic Kondo effect. 



II. MODELS, BARE AND EFFECTIVE 

We consider three semiconducting (single-level) quan- 
tum dots, arranged in a triangular geometry as illustrated 
in Fig 1. Each dot is tunnel-coupled to the others, and 
one of them (dot '2') is also coupled to a metallic lead. 
We focus explicitly on a system tuned to mirror sym- 
metry (see Fig. 1), and study the Anderson-type model 

H = H + H t ri + H hyb . Here H = J^k.a e k a L a k<r refers 
to the non-interacting lead, which is coupled to dot '2' 
via Hhyb = Ek,a v ( a tcr c 2<j + H.c), while H tr % describes 
the isolated TQD with tunnel couplings t,t', 

H trl = ^(e«i + Uh^fiii) 



E [* C 2.K C 



3<T 



H.c. 



(1) 



where hi = h ia = c la c ia 1S the number operator 
for dot i. U is the intradot Coulomb repulsion and 
the level energy of dot i, such that ej = £3 = e for a 




FIG. 1: Schematic illustration of the quantum dot trimer. 
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mirror symmetric system (in which H is invariant under 
a 1 <-> 3 permutation). For convenience we take e, = e 
for all dots (although this is not required, as mentioned 
further below). 

We are interested in the TQD deep in the Af = 3 elec- 
tron Coulomb blockade valley. To this end, noting that 
coupled quantum dot experiments typically correspond 
to t/U <~ 10~ 2 6 , we consider the representative case 
e = —U/2 with U ^> (t,t'). Each dot is then in essence 
singly occupied, with TV = 2 or 4 states much higher 
(~|[/) in energy. The Af = 3 states of the isolated TQD 
comprise two lowest doublets, and a spin quartet (which 
always lies highest in energy). As the tunnel coupling 
t' is increased there is a level crossing of the doublets, 
which are degenerate by symmetry at the point t' = t. 
Projected into the singly-occupied manifold the doublet 
states are 
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with S z = | and cr = ± for spins | / j. Their energy sep- 
aration is E A = E + —E- = J— J' with antifcrromagnetic 
exchange couplings J = 4t 2 /U and J' = 4t' 2 /U, such 
that the levels cross at t' — t (reflecting the magnetic 
frustration inherent at this point). The | — ; S z ) doublet, 
containing triplet configurations of spins '1' and '3', has 
odd parity (— ) under a 1 <-> 3 interchange; while |+; S z ), 
which has singlet-locked spins '1' and '3' and behaves 
in effect as a spin-i carried by dot '2' alone, has even 
parity (+). 

On coupling to the lead the effective model describing 
the system on low energy/temperature scales is obtained 
by a standard Schricffcr- Wolff transformation 12,35 . Pro- 
vided the doublets are not close to degeneracy, only the 
lower such state need be retained in the ground state 
manifold: |-; S z ) for J'« J and |+; S z ) for J' > J. In 
either case a low-energy model of Kondo form arises 



Heff — Jkj S • S(0), 



(3) 



(potential scattering is omitted for clarity), with S(0) 
the conduction band spin density at dot '2', and S 
a spin-i operator representing the appropriate doublet 
(7 = + or — ). The effective Kondo coupling is J^ 7 = 
2(7; +3IS2I7; +\)Jk w hh s 2 the spin of dot '2'; where 
Jk = 8T/(npU) with hybridization T — irV 2 p and p the 
lead density of states. Eqn.(2) thus gives Jk- = —\Jk, 
and Jk+ = +Jk- Hence, for tunnel coupling i' « ! 
(</' <C J), a /erromagnetic spin-i Kondo effect arises 
(Jk- < 0) 25 . Kondo quenching of the lowest doublet 
is in consequence ineffective, and as temperature T — > 
the spin becomes asymptotically free - the stable fixed 
point (FP) is the local moment (LM) FP with a residual 
entropy Si mp (T = 0) = In 2 (k B = l) 33 - The system here 
is the simplest example of a 'singular Fermi liquid' 36 , re- 



flected in the non-analyticity of leading irrelevant correc- 
tions to the fixed point 36 ' 37 . For t' ^> t by contrast, the 
Kondo coupling is antiferromagnetic (Jk+ = +Jk > 0), 
destabilizing the LM fixed point. The strong coupling 
(SC) FP then controls the T — > behaviour, describ- 
ing the familiar Fermi liquid Kondo singlet ground state 
in which the spin is screened by the lead/conduction 
electr ons below the characteristic Kondo scale Tk, with 
T K /VUT ~ exp(-l/pJ K+ ) = exp(-7rl78r) 38 . 

Since the fixed points for the two stable phases are 
distinct, a quantum phase transition must thus occur 
on tuning the tunnel coupling t' through a critical value 
t' c ~ t. We study it below, but first outline the effec- 
tive low-energy model in the vicinity of the transition. 
Here, as the |±; S z ) states are of course near degenerate, 
both doublets must thus be retained in the low-energy 
trimer manifold, and the unity operator for the local 
(dot) Hilbert space is hence: 

i = J2(\+;S z )(+;S z \ + \-,s z )(-,s z \) = + (4) 

The effective low-energy model then obtained by 
Schrieffer- Wolff is readily shown to be 



H. 



= J K is 2 i • S(0) + 5#A(i+ - i 



(5) 



with Jk as above. The final term here refers simply to the 
energy difference between the two doublets. It may be 
written equivalently as EaT z with a pseudospin operator 



t, = i(i + -i_) 



(6) 



thus defined, such that the doublets are each eigenstates 
of it, f z \± ; S z ) = ±i|±; S z ). Considering now the first 

term in eqn.(5), ls 2 l = l + s 2 l + + l_s 2 l_ for the mirror 
symmetric case considered (cross terms vanish by sym- 
metry). Direct evaluation of l±s 2 l± gives l+s 2 l + = S1+ 
(= 1+S) and l_s 2 l_ = — -|l_S, where S is a spin- 
\ operator for the dot Hilbert space (specifically S z — 

E 7 =±, s* \r,s z )s z ( r ,s z \ and5± = E 7 |7;±|><7;t||). 

Hence, using eqns.(4,6) to express 1± = |(1±27^) in 
terms of the pseudospin, the effective low-energy model 
is given from eqn.(5) by 



rrtrans 
U eff 



U K (l+4T z )S-S(Q) + E A f z , (7) 



expressed as desired in terms of the spin S and pseu- 
dospin T z . The term E/\T Z is equivalent to a magnetic 
field acting on the pseudospin, favoring the | — ;S Z ) dou- 
blet for E A > and |+; S z ) for E A < 0; such that cqn.(7) 
reduces, as it should, to one or other of eqn.(3) in the 
limit where the separation \E A \ is sufficiently large that 
only one of the doublets need be retained in the low- 
energy TQD manifold. Finally, note that the absence 
of pseudospin raising/lowering terms T ± in Hlff ns re- 
flects the strict 1 <-» 3 parity in the mirror symmetric 
setup (which cannot be broken by virtual hopping pro- 
cesses between dot '2' and the lead); and means that the 




FIG. 2: Entropy S , imp (T) vs T/T for fixed U = 10 and 
i = 0.1, as a function of t' as the transition (t' c ~ 0.097) 
is approached from either side: for t' = tj. ± XTk with 
A = 10 4 ,10,1,10 _1 ,10 _4 and ~ (curves (a)-(f) respec- 
tively). Solid lines: antiferromagnetically coupled SC phase 
(A > 0); dashed lines: ferromagnetically coupled LM phase 
(A < 0). Inset: close to the transition, the scale Ta vanishes 

~ 2 - 2 

linearly in (£' —t' c ). 

Hilbcrt space of eqn.(7) separates exactly into spin and 
pseudospin sectors, such that only the sign of the effec- 
tive Kondo coupling is correlated to the pseudospin. 



III. RESULTS 

The physical picture is thus clear, and indicates the 
presence of a quantum phase transition as a function of 
t'. We now present NRG results for the TQD Anderson 
model, using a symmetric, constant lead density of states 
p = 1/(2D). The full density matrix extension 30,31 of 
the NRG is employed 34 , together with direct calculation 
of the electron self-energy 39 . Calculations are typically 
performed for an NRG discretization parameter A = 3, 
retaining the lowest N s = 3000 states per iteration. 

As above we choose e — — \U and t/U = 10~ 2 , which 
realistic case 6 corresponds to single occupancy of the dots 
(all calculations give (n;) = 1 for each dot). The low tem- 
perature behaviour is determined by three fixed points: 
those for the two stable phases at T — (SC or LM), 
and a 'transition fixed point' precisely at the transition, 
which at finite-T strongly affects the behaviour of the 
system close to the transition. 

The T-dependence of the entropy 5 l i mp (T) 33 provides 
a clear picture of the relevant fixed points. We show it 
in Fig. 2, for U = U/ttT = 10 and i = t/irT = 0.1 (with 
T/D = 10~ 2 ), for variable t' = t'/nT approaching the 
transition from either side: t' = t' c ± XTk, varying A. 
Here t' c = 0.09715. .(~ t as expected), and the antiferro- 
magnetic Kondo scale T K /T ~ 7 x 10~ 6 38 . Solid lines 
refer to systems in the SC phase (f > t' c ), dashed lines 
to the LM phase (t' < t' c ). In all cases the highest T be- 
haviour is governed by the free orbital fixed point 33 , with 




T/r->- 



FIG. 3: Spin susceptibility Txim P (T) vs T/T for the same 
parameters as Fig. 2. Solid lines: antiferromagnetically cou- 
pled SC phase; dashed lines: ferromagnetically coupled LM 
phase. Close to the transition, for \Ea\ < T < Tk, a 
rximp(r) = | plateau symptomatic of the transition fixed 
point arises, persisting down to T = precisely at the transi- 
tion where Ea = (case (f)). 



all 4 3 states of the TQD thermally accessible and hence 
Simp = ln(64). On the scale T ~ U the dots become 
singly occupied, the entropy thus dropping to ln(8). 

On further lowering T deep in either the LM or SC 
phases (case (a) in Fig. 2), all but the lowest trimer dou- 
blet is projected out and Si mp approaches ln(2), signi- 
fying the LM fixed point. For t' < t' c this remains the 
stable fixed point down to T = 0, while for t' > t' c the 
antiferromagnetic Kondo effect drives the system to the 
SC fixed point below T - T K . Lines (b)-(e) in Fig. 2 
are for systems progressively approaching the transition. 
Here, when T exceeds the energy gap between the dou- 
blets (denoted |£a| and naturally renormalized slightly 
from the isolated TQD limit of |-Ea|)> the pair of dou- 
blets are effectively degenerate 40 and an <Si mp = ln(4) 
plateau is thus reached. The fixed point Hamiltonian 
here is then simply a free conduction band, plus two free 
spins (eqn.(7) with Jk and E\ set to zero). 

For t' within ~ T K of t' c , as in (c)-(e) of Fig. 2, a fur- 
ther decrease in T leads to a clear entropy plateau of 
Simp — ln(3). This is the transition fixed point (TFP): 
for |^ A | < T < T K , the \+;S z ) doublet is screened by 
the antiferromagnetic Kondo effect - even when it is not 
the ground state - while the ferromagnetically coupled 
| — ; S z ) remains a local moment. The TFP thus comprises 
both a free local moment and a Kondo singlet, hence the 
ln(3) entropy. And the TFP Hamiltonian corresponds to 
eqn.(7) with Ea — 0, Jk — ► oo in the T z = +| pseu- 
dospin sector and Jk — > for T z = — ^ ■ The energy level 
spectrum at the TFP thus comprises a set of LM levels 
plus a set of SC levels (as confirmed directly from the 
NRG calculations). 

Finally, on a scale T — Ta ~ |^a|j defined in practice 
by 5i mp (TA) = 0.85 (suitably between In 2 and In 3), Si mp 
crosses over from the TFP value of ln(3) to the T = 



FIG. 4: T = local spectrum ttTD 2 (lo;T = 0) vs lj/F for 
systems with the same parameters as Fig. 2. Panels (A) and 

(B) show all data on logarithmic and linear frequency scales 
respectively, with solid lines again for systems in the SC phase 
(t' > t' c ), and dashed lines for the LM phase (t' < t' c ). Inset 

(C) shows the low-frequency behaviour in the LM phase, fitted 
to a line of the form 7rr£> 2 (^; T = 0) = a/ln 2 (|w|/T ). 



FIG. 5: Conductance G(T)/Go vs T/Tk for the same param- 
eters as Fig. 2. As T — > 0, all systems in the SC phase (solid 
lines) satisfy the unitarity limit G(T — 0)/Go = 1, while 
in the LM phase (dashed lines) G(T = 0)/G = 0. When 
Ta < T < Tk, the conductance is controlled by the transi- 
tion FP, with G(T)/G = § in both phases. 



value appropriate to the stable fixed point (SC or LM). 
As the transition is approached, the scale Ta oc (t n — t 12 ) 
vanishes (Fig. 2 inset) - a natural consequence of the 
doublet level crossing, recalling from above that Ta ~ 
\E A \ = \ J'-J\ = 4\t' 2 -t 2 \/U (with t' c = t for the isolated 
TQD). And since Ta = precisely at the transition, the 
T = entropy at that point is the TFP value ln(3), as in 
case (f) of Fig. 2. 

The behaviour described for Si mp (T) is likewise evident 
in the T-depcndence of Tx imp (T) (with x imp (T) thc to- 
tal impurity/dot contribution to the uniform magnetic 
susceptibility 33 ), as shown in Fig. 3. For T > U gov- 
erned by the free orbital fixed point, Tx imp (T) ~ 3 x | 
as expected 33 (we set g^B = !)• On the scale T <~ U thc 
dots become singly occupied but the spins are essentially 
uncorrelated, so Tx imp (T) ~ 3 x j as expected for three 
free spins 33 . On further decreasing T, the ultimate low- 
temperature behaviour is naturally Tx imp (T) = j for 
the ground state doublet characteristic of the LM phase 
t' < t' c , and Tx imp (T) = for thc quenched SC fixed 
point when t' > t' c . Close to the transition however, for 
|£a| < T < T K , thc TFP is again evident in the per- 
sistence of a Tx imp (T) = g plateau, readily understood 
as the mean ((S z ) 2 } for thc three quasidegenerate states 
arising for \Ea\ < T < Tk as described above; and to 
which value Tx imp (T) tends as T — > 0, precisely at the 
transition (Fig. 3, case (f)). 

Most importantly, the physics above is clearly mani- 
fest in transport properties. The zero-bias conductance 
through the '2' dot is 41 

f°° df(uj) 

G(T)/G = ttT / D 2 (uj;T) dcu, (8) 



with f(ui) the Fermi function and T>2(o>;T) thc local 
single-particle spectrum of dot '2'; such that G(T = 
0)/G = 7rrT> 2 (0;0). Fig. 4 shows thc T = spectrum 
for the '2' dot. Solid (dashed) lines are again for sys- 
tems in the SC (LM) phase. At T = 0, the w = 
spectral density collapses abruptly as thc transition is 
crossed - from the unitarity limit of 7rrZ?2(0;0) = 1 in 
the SC phase (*' > t' c ) to 7rrT> 2 (0; 0) = in thc LM phase 
(f < t' c ). All systems in the SC phase share a common 
Kondo scale, so all solid lines in practice coincide, and 
are found to be characterized by the universal scaling 
form obtained for the single-impurity Anderson model 42 . 
In the LM phase, all T = spectra again coincide. In 
this case however, the low-w behaviour is described by 
ttTD 2 (uj;T = 0) - a/ln 2 (H/T ) (as shown in Fig. 4 in- 
set (C)), as expected for a 'singular Fermi liquid' 37 , in 
which thc slow approach to the fixed point is character- 
ized by marginally irrelevant logarithmic corrections 36,37 . 

The zero-bias conductance as a function of T/Tk is 
itself shown in Fig. 5. At T = 0, G(0)/G o = 1 or in 
the SC or LM phases respectively, as above. Far from 
thc transition (case (a)), G(T) decays steadily with in- 
creasing T in the SC phase, reflecting thermal destruc- 
tion of the coherent Kondo singlet. In the LM phase by 
contrast, G{T)/Gq is not appreciable at any T. Close 
to the transition however, and for Ta < T < Tk, the 
transition FP controls the zero-bias conductance, and 
G(T)/Go = | is seen in both SC and LM phases - mean- 
ing in particular that warming a system in the ferro- 
magnetically coupled LM phase produces a revival of the 
antiferromagnetic Kondo effect. This behaviour is read- 
ily understood by noting that at T = (as in Fig. 4), 
Z? 2 (w;0) = £>2 M (w) f° r t' < t' c , as only excitations from 
the doubly degenerate LM ground states are relevant; 
while for t' > t' c , D 2 (w,0) = D^ \uj) since now only ex- 
citations from thc Kondo singlet arise. At finite-T how- 



■5 



ever, with Ta <C T <C Tr- such that the lowest manifold 
of states comprises both the LM and the Kondo singlet, 
it is easy to show from the Lehmann representation of 
the spectrum that D 2 (cu;T) = \[D^ C {u) + 2D% M (w)] 
(such that ttTD 2 (uj;T) = 1/3 for \u\ < T K ); and hence 
from Eq. 8 that G{T)/Gq = 1/3 - which persists down to 
T = precisely at the transition, where Ta = (Fig. 5, 
case(f)). 

IV. CONCLUSION 

The TQD ring system in the 3-electron Coulomb block- 
ade valley exhibits a rich range of physical behaviour. 
Both antiferromagnetic and ferromagnetic spin-^ Kondo 
physics is accessible on tuning the tunnel coupling t'; 
the two phases being separated by a level crossing quan- 
tum phase transition, reflected in a transition fixed point 
which controls in particular the conductance in the vicin- 
ity of the transition. We also add that while explicit 
results have been given for a mirror symmetric TQD 



with Ci = e (= —£7/2) for all dots, our conclusions 
are robust provided the dots remain in essence singly 
occupied. Varying e (or £7) on dot '2' for example, 
making it inequivalcnt to dots '1' and '3', docs not 
break mirror symmetry and leaves unaltered the be- 
haviour uncovered above. Indeed even breaking mirror 
symmetry, via e.g. distinct tunnel couplings between all 
dots, still results in both ferromagnetically-coupled and 
antifcrromagnctically-couplcd (Kondo quenched) ground 
states separated by a quantum phase transition 43 . The 
robustness of the essential physics suggests that both 
phases should be experimentally accessible in a TQD de- 
vice; as too should the transition between them, provided 
the tunnel couplings can be sufficiently finely tuned. 
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